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(— I \ Abstract 

We demonstrate that so called "infrared divergences" which have been discussed in 
some publications during several years, do not appear within the correct treatment of 
analytical properties of the transition amplitudes, in particular, of the second order pole 
'nJ" ■ structure of the amplitudes describing the n — n transition in nuclei. Explicit calculation 

with the help of the Feynman diagram technique shows that the neutron-antineutron oscil- 
lations are strongly suppressed in the deuteron, as well as in heavier nuclei, in comparison 
f— ^ ■ with the oscillations in vacuum. General advantages and some difficulties of the field 

<^ , theoretical methods applied in nuclear theory are reminded for the particular example of 

the parity violating np — )■ d'y capture amplitude. 
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1 Introduction 

The neutron-antineutron transition induced by the baryon number violating 
interaction [AB = 2) predicted within some variants of grand unified theories 
(GUT) has been discussed in many papers since 1970 [1], see [2] — [6]. Exper- 
imental results of searches for such transition are available, in vacuum (reactor 
experiments [7], and references therein), in nucleus ^^O [8] and in Fe nucleus 
[9], see also the PDG tables. 

During the later time there have been many speculations that the neutron- 
antineutron oscillations in nuclei are not suppressed in comparison with the n—n 
transition in vacuum [10, 11]. The arguments were based on the "true field- 
theoretical approach" to this problem. The result of [10] has been criticized in 
a number of papers [12, 13, 14, 15] which used somewhat different approaches 
(potential, S-matrix, diagram), and general physics arguments. 

However, in view of continuing publications [11] containing same state- 
ment as in [10], it seems to be necessary to analyze this problem just within the 
quantum field theory based approach used in [10, 11]. Our consideration is close 
to the approach of paper [14] where the diagram technique has been applied 
to study neutron-antineutron transition in nuclei, although differs from [14] 
in some details. More recent realistic calculations of the neutron-antineutron 
tansition in nuclei can be found in [16] (the diagram technique motivated con- 
sideration) and in [17] (potential approach). 

In the next section the n — n oscillations in vacuum are considered and no- 
tations used in present paper are introduced. In section 3 we give some general 
arguments based on analytical properties of amplitudes in favour of suppression 
of the n — n transition in nuclei. The simplest example of the deuteron when 
the final result can be obtained in closed form, is considered in details in sec- 
tion 4, where the result of [14] for the case of the deuteron is reproduced. The 
analogy between analytical properties of the amplitude describing the n — n 
transition and the amplitude which corresponds to the nucleus formfactor at 
zero momentum transfer is noted in section 5. The specific difficulties of the 
field-theoretical methods applied to nuclear reactions are recollected for the case 
of the parity violating np — > d^y amplitude in section 6. This concluding section 
contains also some explicit remarks on the approach of [10, 11] and on recent 
E-prints by V.Nazaruk. 



2 The n — n transition in vacuum 

To introduce notations, let us consider first the nn transition in vacuum which 
is described by the baryon number violating interaction (see, e.g. [2, 13, 14]) 
V = finnO'i/2, 0"! bciug the Pauli matrix, finn is the parameter which has the 
dimension of mass, to be predicted by grand unified theories and to be defined 
experimentally ^. As usually, a point-like n — n coupling is assumed here. The 
n — n state is described by the 2-component spinor ^, lower component being 
the starting neutron, the upper one - the appearing antineutron. The evolution 
equation is 

'? = (^» + ^)* (1) 

with Vo = mjM — i^n/'^ in the rest frame of the neutron (mjv is the nucleon 
mass, 7„ - the (anti)neutron normal weak interaction decay width, and we take 
In = Im as it follows from CP-invariance of weak interactions). Eq. (1) has 
solution 



*(t) = exp [-1 {iinntai/2 + Vot)] *o = 






(2) 
Here ^o is the starting wave function, e.g. for the neutron in the initial state 
^0 = (0, 1)"^. In this case we have for an arbitrary time 

^(n, t) = —i sin ™^ exp{—iVot)^ ^(n,t) = cos ™ exp{—iVot), (3) 

which describes oscillation n — n. Evidently, for large enough observation times, 
^obs ^ 1/^^^, the average probabilities to observe neutron and antineutron are 
equal if we neglect the natural decay of the neutron (antineutron): 



W{n) = |*(n)|2 = |*(n)|2 = W{n) = 1/2. (4) 

This case is, however, of academic interest, only, since 7^ ^ Unn ^ It should 
be stressed that in vacuum the neutron goes over into antineutron, also the 



^ There is relation /i„,-i = 2Sm with the pararaeter Sm introduced in [2]. The neutron- antineutron oscillation 
time in vacuum is Tnn = l/Sm = 2/ finn, see also [17] and references in this paper. 

^It is a matter of simple algebra to calculate the integrals over time of the probabilities |^(n,t)p and 
l*(n,i)p:. 
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for neutron as initial state and for arbitrary, but different from zero 7„. The difference between both quantities 
is obvious, and disappears when 7„ — > 0. 



discrete localized in space state, which can go over again to the neutron, so the 
oscillation neutron to antineutron and back takes place. 

Since the parameter finn is small, the expansion of sin and cos can be 
made in Eq. (3) at not too large times. In this case the average (over the time 
^obs ^ 1/^^^) change of the probability of appearance of antineutron in vacuum 
is (for the sake of brevity we do not take into account the (anti)neutron natural 
instability which has obvious consequences) 



.2 j-obs 



which has, obviously, dimension of the width F. So, in vacuum the transition 
n — > n is suppressed if the observation time is small, t"^^ <C l/l^nn- From 
existing data obtained with free neutrons from reactor the oscillation time is 
greater than 0.8610^sec ~ 2.7 years [7], therefore, 

linn < 1.510-23 eV, (6) 

very small quantity. 

Recalculation of the quantity finn or r„^ from existing data on nuclei 
stability [8, 9] is somewhat model dependent, and different authors obtained 
somewhat different results, within about 1 order of magnitude, see e.g. dis- 
cussion in [14, 16, 17]. Most recent results for finn obtained from the nuclear 
stability data are close to (6) [16, 17], see also the next section. 

3 Analyticity based arguments for the suppression of the 
n — n transition in a nucleus 

In the case of nuclei the n — n line with the transition amplitude finn is the 
element of any amplitude describing the nucleus decay A ^ (^4 — 2) + mesons, 
where (^—2) denotes a nucleus or some system of baryons with baryonic number 
A — 2, see Fig. 1. The decay probability is therefore proportional to fi^f^, and 
we can write by dimension arguments 

r(A ^(A-2) + mesons) - ^, (7) 

mo 

where mo is some energy (mass) scale. For the result of [10, 11] to be correct, 
the mass mo should be very small, mo ~ finn ~ lO^^^ eV, but we shall argue 



that iJio is of the order of normal hadronic or nuclear scale, mo ~ TUhadr ~ 
(10 — 100) MeV. We can obtain the same result from the above vacuum formula 
(5), if we take the observation time t"^^ ~ l/rrihadr- 
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Figure 1: The Feynman diagram describing the n — n oscillation in a nucleus A with subsequent 
annihilation of antineutron to mesons. The final state has the baryon number A — 2. 



Indeed, the matrix element of any Feynman diagram containing such tran- 



sition 



finn{A-Z) V{A-n,{A-l)) 



-2m{A 



T{A ^ {A-2) + mesons) ~ 
f{n+{A-l)^ {A 



2) + mesons) 



(E„ - E» + isy- 



zfJyiliE. ^ K» 



dEr. 



n) ' 



dE„ ~ 



(8) 



according to the Cauchy theorem known from the theory of functions of complex 
variable. E^ is the neutron (antineutron) energy - integration variable, E^ is the 
(anti)neutron on-mass-shell energy E"^ ~ m]s[+jf/2mN. The energy-momentum 
conservation should be taken into account for the vertex V{A -^ n + {A — 1)) 
which includes the propagator of the (^4 — 1) system, and for the annihilation 
amplitude T. The case of the deuteron considered below is quite transparent 
and illustrative. 

The amplitude T which describes the annihilation of the antineutron, and 
the vertex function V are of normal hadronic or nuclear scale and cannot, in 
principle, contain a very small factors in denominator (or very large factors, of 
the order of 10^^, in the numerator). By this reason we come to the above Eq. 
(7), and the resulting decay width of the nucleus is very small, 

(9) 



30, 



r(^ ^ (^ - 2) + mesons) < lO^''^, 



nnj 



at least 30 orders of magnitude smaller than the inverse time of neutron- 
antineutron oscillation in vacuum ^nn- From Eq. (7) or (9) we obtain 



^^nn ~ \J^{'^ ~^ (^ ~ 2) + mesons)mQ, (10) 

and when one tries to get the restriction on finn from the data on nuclei stability 
[8, 9] the result is close to that from the vacuum experiment [7], somewhat 
smaller, within one order of magnitude [5, 6, 14]. The result of [16] based on 
the intuitive physical picture of n — n transition in medium, differs from that 
of [14] for nuclei ^^O and Fe, and the authors [16] come to the conclusion, that 
experiments with free neutrons from reactor could provide stronger restriction 
on the neutron-antineutron transition parameter than experiments on stability 
of nuclear matter ^. 

According to [10, 11] the probability of the nucleus decay is proportional 
to W{t"^^) ~ f4ini't"^^) (the process proceeds similar to the vacuum case). 



where t°^^ is the large observation time, of the order of ~ 1 year or greater. 
By this reason the extracted value of finn is smaller than that given by Eq. 
(10), by about 15 orders of magnitude. Technical reason for strange result 
obtained in [10, 11] is the wrong interpretation of the second order pole structure 
of any amplitude containing the n — n transition. Instead of using the well 
developed Feynman diagram technique, the author [10, 11] tries to construct 
the space-time picture of the process by analogy with the vacuum case, which 
is misleading. Further discussion of papers [10, 11] and recent E-prints of this 
author can be found in concluding section 6. 

4 The case of the deuteron 

We continue our consideration with the case of the deuteron which is quite 
simple and instructive, and can be treated using the standard diagram technique 
^. The point is that in this case there is no final state containing antineutron 
— it could be only the pn state, by the charge conservation. But this state is 
forbidden by energy conservation, since the deuteron mass is smaller than the 



•^There is, in fact some kind of competition between both methods, and final resuh wiU depend on the 
progress to be reached in both branches of experiments — with free neutrons and with neutrons bound in 
nuclei. Friedman and Gal [17] obtained the restriction Tnn > 3.310^sec from the latest datum on ^^O stability 
and using the potential approach. Experiments with ultracold neutrons in a trap have been proposed and 
discussed in [3, 18], but not performed till now. 

^It has been considered in fact in [19] within the reasonable framework of the diagram technique. However, 
the author has drawn later wrong conclusions from this consideration. 



sum of masses of the proton and antineutron. Therefore, if the n — n transition 
took place within the deuteron, the final state could be only some amount of 
mesons. The amplitude of the process is described by the diagrams of the type 
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Figure 2: The Feynman diagram describing n — n oscillation in the deuteron with subsequent 
annihilation of antineutron and proton to mesons. 



shown in Fig. 2 and is equal to 

T{d -^ mesons) = igdnpruNldnn 



T{np 



mesons. 



d'^p 



{p^ - m%)[{d - pf - mlY (27r)4" 
The constant gdnp is normalized by the condition [20, 21, 22] 



ir 



9dnp _ ^ _ €rf 

IGtt m^r \ niN 



0.049, 



(12) 



which follows, e.g. from the deuteron charge formfactor normalization Ffi{t = 
0) = 1, see the next section, k = s/niN^d, Q — 2.22 MeV being the binding 
energy of the deuteron. For the vertex d ^ npwe are writing 2mNgdnp to ensure 
the correct dimension of the whole amplitude, see also the next section 

The integration over internal 4-momentum d'^p in (11) can be made easily 
taking into account the nearest singularities in the energy po = E", in the non- 
relativistic approximation for nucleons. As we shall see right now, the integral 
over d^p converges at small j9 ~ k which corresponds to large distances, r ~ 1/k. 
By this reason the annihilation amplitude can be taken out of the integration 
in some average point, and we obtain the approximate equality 



T{d — > mesons) = gdnpTnNfJ'nnIdNNT{np — > mesons) 



(11a) 



with 

IdNN = 



{27TyJ {p^ - m%)[{d - py - m%]^ 

d'^p 



(27r)4(2m)3 J [p^ - jun - p^/{2mN) + i6] [nid - tun - Po - p^/{2mN) - iS]"^ 

d^p 1 



This integral converges at small \p\ ~ k, more details can be found in the next 
section: the integral IdNN enters also the deuteron charge formfactor at zero 
momentum transfer. The decay width (probability) is, by standard technique, 

T(d — > mesons) ~ f^nn9dnp-^dNN''^N / \T{fip — > 'mesons)\ d^{mesons), (14) 

^{mesons) is the final states phase space. Our final result for the width of the 
deuteron decay into mesons is 



167r/t'"" ^'"' ^■-/.^;o^0 g^^ 



Per. 



where Pc.m. is the (anti)nucleon momentum in the center of mass system. This 
result is very close to that obtained by L.Kondratyuk (Eq. (17) in [14]) in 
somewhat different way, using the induced np wave function ^. 

The annihilation cross section of the antineutron with velocity vq on the 
proton at rest equals 

(T{np — > mesons) = 2 — / \T{np -^ mesons)\'^ d^ (mesons) . (16) 



4mj^vo 

According to PDG at small vq, roughly, vga^^ ~ (50 — 55)m6 ~ (130 — 

UQ)GeV^\ So, we obtain from (15) /inn < 2.510-24e"l/, or Tnn > 5T0Sec if 
we take optimistically the same restriction for the deuteron stability as it was 
obtained for the Fe nucleus, r^ ~ rpe > G.blO^^yr [9]. Our result (15) is valid 
up to numerical factor of the order ~ 1, since we did not consider explicitly 
the spin dependence of the annihilation cross section and the spin structure of 



^The result Eq. (17) in [14] can be rewritten in our notations as 

2 

which differs from our result by some numerical factor, close to 1 and not essential for our conclusions. 



the incident nucleus. Same holds in fact for the results obtained in preceeding 
papers, see e.g. [2, 14]. 

Additional suppression factor in comparison with the case of a free neu- 
tron is of the order of 

in agreement with our former rough estimate (9), and disappears, indeed, when 
the binding energy becomes zero ^ . The binding energy of the deuteron should 
be very small, to provide the value n ~ ^nn', to avoid such suppression. At such 
vanishing binding energy the nucleons inside the deuteron are mostly outside 
of the range of nuclear forces, similar to the vacuum case. 

Results similar to (15) can be obtained for heavier nuclei, see [5, 6, 14, 16, 
17]. The physical reason of such suppression is quite transparent and has been 
discussed in the literature long ago (see e.g. [2, 13, 15]): it is the localization 
of the neutron inside the nucleus, whereas no localization takes place in the 
vacuum case. In the case of the deuteron or heavier nucleus the annihilation 
of antineutron takes place, and final state is some continuum state containing 
mesons. By this reason the transition of the final state back to the incident 
nucleus is not possible in principle, and there cannot be oscillation of the type, 
e.g. d — > mesons — > d. This is important difference from the case of the free 
neutron. 

5 The deuteron charge formfactor 

As we noted previously, the presence of the second order pole in intermediate 
energy variable is characteristic for the processes with the neutron - antineutron 
transition, but it is in fact not a new peculiarity, it takes place also for the 
case of the nucleus formfactor with zero momentum transfer, FA{q = 0). Let 
us consider as an example the deuteron charge formfactor. In the zero range 
approximation it can be written as 



i{2 mgdnpf r d^ 

(27r)4 J (^2 _ ^2^) [(^^ -p)^ - m%] [{d-p + qY - m\] 



^d[q) - — 77TZU — / 7Z9 — _2 ATM — :^2 — _2ir/j — :r—K9. — iryr- i^^) 



Behind the zero range approximation g^np should be considered as a function of 
the relative n — p momentum, not as a constant. For q = second order pole 



®There is no final formula, for Td^mesons in [19] to be compared with our result (14), (15). Numerically, 
however, the result of [19] is in rough agreement with our and [14] estimates. 

9 



appears, and we come to the expression for F{q = 0) containing the integral 
IdNN introduced above in Eq. (13): 



i2 



Fd{0) = {2mN9dnpYhNN. (18) 




Qdnp Qdnp 

Figure 3: The Feynman diagram describing the deuteron charge formfactor. 

In the nonrelativistic approximation, when only the nearest in energy 
E = pq singularities are taken into account, the integral over the energy has the 
structure 

'""' ^ J (E-a + rS)(E-b-rSf = (^^6)^' ^^^^ 

a = rriN + ff/2mN. b = m^ — ijin — ^/2m7v, a — 6 = e^ + ff/rriN, and can be 
calculated using the lower contour which includes the po[e at E" = a — id, or the 
upper contour, including the second order pole at E = b + id, with the help of 
formulas known from the theory of functions of complex variables. After this 
we obtain 



Since Fd{0) = 1, this leads to the above mentioned normalization condition 
fi'lp/(167r) = ^Jcd/mN ^. 

This relation between the constant gdnp and the binding energy of the 
weakly bound system (deuteron in our case) is known for a long time [20, 21, 22]. 
It was obtained in [20, 21, 22] using different methods, dispersion relation, for 

^As it is known from the nonrelativistic diagram technique, the wave function of the deuteron in momentum 
representation is ^(i(p) = 5dnp/ [47r'^/^(K^ +p^)], therefore, the normahzation of the charge formfactor Frf(O) ~ 1 
foUows from the normahzation of the deuteron wave function, which is also well known from quantum mechanics. 

10 



example. We shall demonstrate here for completeness, following to Landau 
[22], that relation (12) between the constant gdnp and binding energy appears 
from the consideration of the pole contribution to the two-particle scattering 
amplitude, the n]9-scattering in our case, see Fig. 4. 




Vdnp ^ 9dnp 




Figure 4: The Feynman diagram corresponding to the deuteron pole in the np scattering 
amphtude 



The contribution of the pole diagram (Fig. 4) to the relativistic invariant 
scattering amplitude due to the presence of the bound state (the deuteron in 
our case) equals 

(21) 



rppole ^ {^rUNgdnp 



np^np 2 ' 

s — m^ 

where the Mandelstam variable s = {pn + PpY- At the threshold, s = 4m^, we 
have 

TZ%n,{s = 4mi,) = !^. (22) 

since at the threshold s — m\ = AwiN^d = 4/v^ and we assume for simplicity that 
both the proton and neutron masses are equal to m^y. 

Now we should compare this result with the known quantum-mechanical 
expression for the scattering amplitude in the zero range approximation 

k being the value of the nucleon 3-momentum in the center of mass frame. Using 
the known relation between the relativistic invariant and quantum-mechanical 
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scattering amplitudes, T{s) = 87^^/sf{k), at the threshold {k = 0) we obtain 

(24) 



^N9dnp _ IfiTrmjv 



and relation (12) follows from (24) immediately. 

If the infrared divergence discussed in [10, 11] took place for the process 
of n — n transition in nucleus, it would take place also for the nucleus formfactor 
at zero momentum transfer. But it is well known not to be the case, as we also 
illustrated in this section. 

6 Concluding discussion and remarks 

The field-theoretical description of nuclear reactions and processes is potentially 
useful, it allows to study some effects which is not possible, in principle, to study 
in other way, e.g. relativistic corrections to different observables. One should be, 
however, very careful to treat adequately analytical properties of contributing 
amplitudes. 

In the case of the parity violating amplitude of the radiative capture 
of the low energy neutrons by protons relativistic contributions change the 
nonrelativistic weak interaction isospin selection rules for the parity violating 
observables: photon circular polarization (neutrons unpolarized) and photon 
asymmetry in the capture of polarized neutrons. This has been a motivation 
to study such relativistic contributions to parity violating observables in the 
np ^ d'y - reaction [23]. In this case it was necessary to take into account 
contributions of all singularities (poles) of the amplitude in the complex energy 
plane of the virtual nucleon, not only contributions of the nearest poles in the 
energy variable, as it is made usually in the nonrelativistic calculations. Besides, 
and it is the spesifics of the processes with photon emission, the contact terms 
should be reconstructed to ensure the gauge invariance of the whole amplitude 
of the photon radiation [23] . The nonrelativistic diagram technique developed 
up to that time turned out to be misleading for the case of physics problem 
considered in [23]. The cancellation between contributions of different poles 
has not been noted in first publications on this subject [24]. As a result of 
this cancellation the relativistic contributions to the observables turned out to 
be not greater than nonrelativistic values, in spite of the change of the isospin 
selection rules. 

12 



This particular example is only one of many possible illustrations of the 
difficulties of field-theoretical methods applied to various nuclear physics prob- 
lems. It is often not so easily and straighforwardly to resolve appearing con- 
tradictions with widely known methods and results, as it happened in the case 
of papers [10, 11]. The author of [10, 11] tries to reconstruct the space-time 
picture of the process, but the correspondence of this picture to the well jus- 
tified amplitude, as it appears from the Feynman diagrams, is questionable. 
The infrared divergence discussed in [10, 11] is an artefact of this inadequate 
space-time picture of the whole process of n — n transition with subsequent 
antineutron annihilation. A reasonable and logically consistent way would be 
to rewrite the amplitude which corresponds to Feynman diagrams with second 
order pole in energy-momentum variables (Fig.l and 2) in space-time variables 
which will provide the correct space-time picture of the process, instead of writ- 
ing ad hoc the amplitude in space-time variables similar to that of the process 
in vacuum. Another quite unrealistic consequence of this space-time picture 
[10, 11] is the nonexponential law of the nucleus decay. There is no "new limit 
on neutron - antineutron transition" [10]; instead, one should treat correctly 
singularities of the transition amplitudes in the complex energy plane. 

In his comment [25] Nazaruk makes the statement: "For the propagator 
in the loop the infrared divergence (for nfi transition, nucleus formfactor and 
so on) cannot be in principle. In order to obtain the infrared divergence the 
neutron line entering the nfi transition vertex should be the wave function." It 
means that the author of [25] agrees that within the Feynman diagram technique 
the "infrared divergence" does not appear, but new rules seem to be proposed in 
[10, 11, 25] instead of well known Feynman rules. These "new rules" should be, 
at least, clearly formulated, and, second, these rules should allow to reproduce 
all well known results of nuclear theory. In [25] there is also the statement 
concerning section 3 of [26] (preliminary version of present paper): "The main 
statement of this section is completely wrong" . However, there are neither 
proof, nor scientific arguments that our results are invalid. 

In recent E-prints [27] the upper bound for the free-space n — n oscillation 
time is extracted from existing nuclear data to be ~ 10^^ years. This is the 
result of the so called "model with bare propagator" and repeats the previous 
statements of [10, 11]. We have just shown in present paper that calculations 
made in [10, 11] according to this "model with bare propagator" are wrong. At 
the same time, the questions put in [26] and here, are not answered in [27]. 

13 
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